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ABSTRACT 

We derive constraints on the relic abundance of a generic particle of mass 
~ 1 — 10 14 TeV which decays into neutrinos at cosmological epochs, using data 
from the Frejus and 1MB nucleon decay detectors and the Fly's Eye air shower 
array. The lifetime of such unstable particles which may constitute the dark matter 
today is bounded to be greater than ~ 10 14 — 10 18 yr, depending on the mass. 
For lifetimes shorter than the age of the universe, neutrino energy losses due to 
scattering and the expansion redshift become important and set limits to the ability 
of neutrino observatories to probe the early universe. 



1. Introduction 



Upper limits on the flux of high energy cosmic neutrinos obtained from nucleon 
decay experiments and cosmic ray observatories constrain the relic cosmological 
abundance of heavy unstable particles which decay into neutrinos. Given the 
energy spectrum of the decay neutrinos and the decay branching ratio, upper 
bounds can be obtained on the primordial abundance of the particle as a function 
of its lifetime. Earlier attempts to set such bounds [1,2] were made before any 
experimental data were available. We present here the bounds imposed by the non- 
observation of extraterrestrial high energy neutrinos in the Frejus and 1MB nucleon 
decay detectors and the Fly's Eye air shower array. We improve on previous work 
by taking into account the experimental energy thresholds, the neutrino opacity of 
the early universe, neutrino absorption in the Earth and the appropriate neutrino 
interaction cross sections at high energies. 

In section 2, we discuss the cosmological absorption of high energy neutrinos, 
and in section 3 calculate the spectrum of neutrinos generated by heavy particle 
decay. A discussion of the expected signals is given in section 4 and the constraints 
provided by present observations are presented in section 5. Our conclusions follow 
in section 6. 

2. Cosmological neutrino absorption 

High energy neutrinos can be absorbed in interactions with the relic ther- 
mal neutrino background and with nucleons in the early universe. The dominant 
processes are the annihilation of a high energy neutrino (or antineutrino) with a 
background antineutrino (or neutrino) and its inelastic scattering off a nucleon. 
We obtain below an analytical formula for the absorption redshift z a (E e ) at which 
the neutrino opacity of the universe s v is unity for a neutrino emitted with energy 
E e ; less than a fraction 1/e of the neutrinos emitted at redshifts larger than z a (E e ) 
propagate to the present epoch. 
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Consider a neutrino emitted with energy E e at time t e corresponding to redshift 
z e . The cosmological neutrino opacity s u (t e , E e ) is the mean number of scatterings 
undergone by the neutrino in which it could have been absorbed, given by 



v(^e, E e ) — J 



<-., 

to 



where to ~ 0.65 x 10 10 yr (QqIi 2 )^/ 2 is the present age of the universe and r u (t, E v ) 
is the mean free time between collisions at time t (and redshift z) for a neutrino 
of energy E v = E e (l + z)/(l + z e ). Here Qq is the present mass density of the 
universe in units of the critical density p c ~ 1.9 x 1CT 29 /i 2 g cm -3 , where h is the 
Hubble constant in units of 100 km s" 1 Mpc -1 (0.4 &h£, 1). Taking into account 
the two absorption processes mentioned above: 



111 . . 

(2.2) 



T v (t,E v ) r V y r^ N 



i 



where the first term refers to neutrino-antineutrino annihilation and the second to 
neutrino-nucleon scattering. 

To obtain t^, we must average over the thermal energy distribution of the 
relic background (anti)neutrinos. Consider a decay neutrino and a background 
antineutrino; the same formulae apply to a decay antineutrino and a background 
neutrino. Indicating by 9^ the angle between the two colliding particles in the 
cosmic frame, we have 



— = ((1 - cos p ^)a V y)n^ , (2.3) 



where (1 — cos 6^) is the vv relative velocity (in units of c), a V y is the total vv 
annihilation cross section and rip is the background antineutrino number density 
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at temperature Tp, given by 



n- - ^T 3 

V ~ A 7T 2 V ' 



(2.4) 



The angular brackets indicate an average over the antineutrino energy distribution, 



i 



2vr 2 q e,/t, + x • 



(2.5) 



We consider only annihilations into charged fermion pairs, vv — > //. In this 
case, o vv = J2f a vv^fji w ith the sum running over quarks and charged leptons. It 
is a good approximation for our purposes to consider massless fermions to estimate 
the annihilation cross section and simply add a new fermion channel whenever 
E v Ty > mj. In this case 



An 



cc 



(2.6) 



with s = 2EyE u (l — cos9yy). The first term in the square brackets is due to neutral 
currents. The Z boson pole factor is defined as 



P 7 ,(s) 



(s - M|)2 + M\Fl ■ 



(2.7) 



The second term takes into account the charge current contribution to the processes 



v e v e ->• e+e , v^v^ -> pL+pL , v T v T -> r+r : 



A w (s) 



with 



Mir 



2s 3 sin #\y 



. . s as 2 , .f s 

(2.8) 



(2.9) 



(s - M|)2 + J?Af| • 

The coefficients iV e ff are the effective numbers of annihilation channels. For neutral 
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currents, this is calculated as 

^eff = J2 9 ( E " T " ~ m /) hfi 1 - 8t 3f1f sin2 d W + 8qj sin 4 9 W ) , (2.10) 
/ 

while for charged currents, the coefficient 

Np s c = 9 {E v Ty - mf) f sin 2 9 W (2.11) 

is non-zero only if the charged lepton I is in the same family as the annihilating 
neutrino. Above, rif is the number of colours (1 for leptons, 3 for quarks), and 
t^f and qf are the third component of the weak isospin and the electric charge of 
the fermion in units of the positron charge respectively. We take the electroweak 
mixing angle to be given by sin 2 = 0.23. Inserting eq. (2.6) into eq. (2.3), we 
obtain 

^T = 3( (1_C0SW [ N eS C sPz(s) + N^ s c sA w (s)\ ^n„, (2.12) 

where s = 2E v Ev(l — cos#^) is understood. 

The thermal average has a peak at E V T„ ~ M|/4, corresponding to the Z pole. 
All neutrinos emitted with energy E e ^ M|/4TV = 1.26 x 10 13 TeV/(l + z e ) are 
absorbed. For E e & M^/AT^, the factors Pz(s) and A^(s) in eq. (2.6) can be set 
equal to unity and is easily evaluated as 

— = <T N eS E u p^ = p^a G N e ^ l + Z ^ E e . (2.13) 

Here, iV e fj = N^p + N^p is the effective total number of annihilation channels 
(which varies between 0.33 and 10.1), <jo is defined as 

<7q = —Gl = 1.12 x 10~ 32 cm 2 TeV" 2 , (2.14) 

and 

p Po = 5.96 x 10- 14 TeVcm~ 3 (2.15) 
is the present antineutrino energy density (per species), taking the present photon 
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temperature to be 2.74 K [3]. Thus, the annihilation mean free time is 



r vV = 1.58 x 10 27 yr N^(l + *)-5(i + Ze ) . (2.16) 



Next we consider neutrino-nucleon scattering. The thermal motion of the non- 
relativistic nucleons can be neglected, and we have 



1 cr vN (1 + z) 4 



TvN E v 1 + Z e 

where un is the nucleon number density at redshift z and ovn is the neutrino 
nucleon scattering cross section at neutrino energy E v . The present nucleon mean 
density is in the range n^ ~ (0.25 — 1.5) x 10~ 7 cm~ 3 according to Big Bang 
nucleosynthesis calculations; this reflects the observational uncertainty in primodial 
^He mass fraction, which is taken to be ~ 0.21 — 0.24 [4]. 

For E u &lTeV, the ratio o v ^jE v is constant and equal to 0.67 x 10 _35 cm 2 
TeV -1 for neutrinos and to 0.34 x 10~ 35 cm 2 TeV -1 for antineutrinos [5-6]. The 
neutrino scattering mean free time at these energies is 

r, N ^ 10 24 yr(l + ^)- 4 (l + z e )(^ \ (2.18) 

Comparing with the annihilation mean free time (2.16), we see that inelastic scat- 
tering upon nucleons dominates only at redshifts 1 + z&lQpNZfi . However, now 
r v ~ t^n *5 10 15 yr ^> to, i.e. the universe has already become transparent to neu- 
trinos. At higher neutrino energies, a v ^jE v decreases and neutrino-nucleon scat- 
tering is even less important, becoming negligible at all redshifts for E v ^ 10 6 TeV. 
Thus inclusion of z/N scattering affects z & (E e ) only marginally. For simplicity of 
presentation, we do not therefore write it explicitly in the formulae below, although 
we have included it in the numerical calculations. 
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The last ingredient necessary to compute the neutrino opacity is the relation- 
ship between the age of the universe and the redshift: 



(o(1+Jc,) 1/2 (1 + j)" 2 , for 

where l + £ e q = 2.25 x 10 4 f2o/i 2 is the redshift at which the energy density of matter 
(with present density parameter Qq) begins to dominate over that of radiation. 

The absorption redshift obtained from integration of eq. (2.1), using eqs. (2.13), 
(2.17) and (2.19), is shown as the diagonal full line in the E e — t e (or E e — z e ) plane 
in figure 1 taking QqIi 2 = 1. The dot-dashed line separates the two regions where 
annihilation and scattering absorption dominate. The region of interest extends 
from the present epoch (z = 0), through the epoch of matter-radiation equality 
(z~2x 10 4 ), up to the epoch of light neutrino decoupling [z ~ 10 10 ), which are 
all indicated. The location of the Z boson pole is also shown as a diagonal dashed 
line. 

Approximate expressions for the absorption redshift z a (E e ) can be obtained for 
1 < z e < z eq and z e ^> z eq . In these cases, the result of the integration simplifies 
to 



3.5 x l(T 17 W 2 r 1/2 (l + ^c) 5/2 (E c /TeV) , for 1 < z e < z eq , 
0.81 x 1(T 14 (1 + z e ) 2 (£ e /TeV) , for z c > z eq . 

The absorption redshift z a (E e ) is then obtained by setting s v — 1: 

. 3.8 x 10 6 (J?o^ 2 ) 1/5 (£ e /TeV)~ 2/5 , £ e £5.2 x 10 5 TeV(Q h 2 )- 2 , 
1 + z & (E e ) = 4 

1.1 x 10 7 (Ee/TeV) 7 , E e £ 5.2 x 10 5 TeV(A^ 2 ) -2 - 

(2.21) 
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3. Neutrino spectrum 



We now determine the present energy spectrum of neutrinos originating from 
the decay of an unstable heavy particle x with decay lifetime r x . The number of 
neutrinos of type V{ {i>i = u e , P e , u^, 9^, . . .) produced at time t, per unit comoving 
volume and unit time, is 

B v Y x (t) B U Y X f t\ B V Y X ( t-t Q \ 
7e (0 = — = — — - expl -- = — — - expf — , (3.1) 

>X >X \ IX/ >X \ I X / 

where B Vi is the number of neutrinos of type V{ produced per decaying x particle, 
Y x (t) = n x (t) /n-y(t) is the x particle number density in ratio to the thermal photon 
density n^(t) (= 412.7(1 + z) 3 cm -3 ), Y Xp is its primordial value* and Y Xo is its 
value today. 

The number of neutrinos absorbed in the same volume, 7 a (£), is proportional 
to the comoving density of decay-generated neutrinos Y Vi {t) = n Vi {t)/n^{t) and is 
given by 

7a(t) = , (3.2) 

where r Vi {t,E G ) is the neutrino absorption mean free time (see section 2). 
The evolution of the comoving neutrino density Y Vi (t) is governed by 
dY Vi (t) 



dt 



= 7e (t)-7a(t) = 

BuXx v ( t\ Y Vi (t) 
exp — 



(3.3) 



Tx V T xJ T Vi (t,E e ) 

with the following solution at the present epoch to : 

to 

Y Vln = B V Y X \ exp 



^Y Xp J - 
o 



T x 



dt 



e 
T x 



(3.4) 



Now differentiating with respect to the present neutrino energy E UiQ = E e (l + 



★ This is conveniently measured at the earliest epoch following which the comoving photon 
number is conserved, say at T ~ 0.01m c , corresponding to t ~ yr; this is negligible 

compared to all other time-scales relevant here. 
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z e ) 1 , one obtains the present neutrino flux 

t, 



E n tt^- = <t> l0 B v Y Xv k ^ exp 



9(E e -E Uio ), (3.5) 



where 



. . (E,„X J 2, for < £ e (l + 2eq )-\ 

t - = t0 Uv ' K= \\, tor E^> m+ ^-\ (36) 

and 4> l0 = n 7o /47r = 0.98x 10 12 cm~ 2 s _1 sr _1 is the present background photon flux 
per unit solid angle. The decay neutrino flux is shown in figure 2 for r x = 10~ 5 to 
and 3E = 10 5 TeV, 10 7 TeV, 10 9 TeV. Notice that the present neutrino energy 
E v is redshifted from E e . The dotted lines indicate what the flux would have 
been without cosmological neutrino absorption. These three curves are simple 
translations of each other, since the differential flux (3.5) with s Vi = depends 
only on the ratio E ViQ /E B . 

Approximating the effect of the cosmological neutrino absorption with e _Sl/ * ~ 
9(t e — t a ), where t a < to corresponds to the absorption redshift z a (E e ) at which 
the neutrino opacity is unity, eq. (3.5) can be easily integrated to obtain the total 
neutrino flux today, 

<t> Vio ~ 7oJ B^ p (e-^ - e-*°/ T -). (3.7) 

For r x <^to — t a , this reduces to 

4> ViQ ^4> 10 B Vi Y Xp e-^ T % (3.8) 

while for r x ^> to — 1&, it becomes 

^-^B^Yj^l. (3.9) 

The neutrino flux is exponentially suppressed for t x <C t a and reaches a maximum 
of 4>v lQ ~ B Vi Y Xp x 10 12 cm -2 s _1 sr _1 for r x ~ to — t&- This flux is potentially 



9 



enormous compared with present bounds on the diffuse extragalactic high energy 
neutrino flux (~ f CP 6 cm~ 2 s _1 sr _1 for E v ^ 1 TeV and ~ 10~ 16 cm -2 s _1 sr _1 
for E v k,10 7 TeV) which may be inferred from data obtained with underground de- 
tectors and cosmic ray observatories (cf. section 4). Hence very restrictive bounds 
may be obtained on the abundance of the hypothetical decaying particle as demon- 
strated below. 

4. Expected signals 

At present, the best means to detect a diffuse background of high energy neutri- 
nos is through the production of an energetic charged lepton in the collision of such 
a neutrino with a nucleon. We consider three possible types of signal according to 
where the interaction occurs. An event is called contained when the interaction 
occurs inside an underground detector, such as Frejus, 1MB and Kamiokande. A 
flux of through-going muons is registered when interactions occur in the material 
surrounding the detector (rock in the underground experiments and water in the 
forthcoming DUMAND, GRANDE and Lake Baikal experiments). Finally, if the 
interaction occurs in the atmosphere an extensive air shower (EAS) is generated, 
which can be detected by cosmic ray observatories such as Fly's Eye and CASA. 

We consider the following experimental constraints on the total neutrino flux: 

(1) the rate of contained events in the Frejus detector, with electron and/or 
muon energies greater than 3 GeV, does not exceed 17.7 kton _1 yr _1 [12];* 

(2) the rate of contained events with energies between 100 MeV and 2.5 GeV 
in the IMB-3 detector is limited by 111.5 kton _1 y r_1 [13]* 

★ We consider the 11 electron and 14 muon charged current events over 3 GeV observed in 
1.56 kton yr (see fig. 3 of ref. [12]), and apply the quoted identification efficiencies of 85% 
and 95% respectively. 

t From the total number of 422 contained events in 3.4 kton yr, we exclude the 43 events 
below 100 MeV (see fig. 2 of ref. [13]) where the track reconstruction and identification 
efficiencies are low. For comparison the IMB-1 detector recorded 401 contained events in 
3.77 kton yr [14]. 
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In fact, the observed contained events are well accounted for by the expected 
neutrino flux from cosmic ray interactions in the atmosphere [15], within the un- 
certainty of ~ 25% in these computations. Hence the bound on contained events 
of non- atmospheric origin can, in principle, be improved by up to a factor of ~ 10 
and the limits to be derived strengthened proportionally. 

We also consider the following constraints on any extraterrestrial neutrino flux: 

(3) the flux of upward-going muons (from directions with zenith angle larger 
than 98°) with energy greater than 2 GeV registered by the IMB-1 detector is less 
than 2.65 x 1CT 13 cm~ 2 s _1 at the 90% confidence level, after subtraction of the 
expected atmospheric component [16]; * 

(4) the Fly's Eye array has set upper limits on the rate of neutrino-induced 
EAS's of 10~ 45 s- 1 sr-\ 3.8x lO" 46 s" 1 sr" 1 , 10~ 46 s" 1 sr" 1 and 3.8x 10~ 47 s" 1 sr" 1 , 
all at the 90% confidence level, for neutrino energies higher than 10 5 TeV, 10 6 TeV, 
10 7 TeV and 10 8 TeV respectively [18]. 

For the isotropic neutrino flux (3.5), the rate of contained events per unit 
detector mass, R c , the upward-going muon flux, 0^, and the rate of EAS's per unit 
solid angle, J, can all be written in the form 



where the signal S is R c , 0^ or J, and the sum is over neutrino types (z/j = v e , 
P e , z/ M , 9^, . . .). The effective aperture Qi(E Ui ), which takes account of neutrino 
absorption by the Earth (if any), and the transfer functions Pi(E Vi ) depend on the 
experimental data set considered. (We have dropped the subscript referring to 
the present neutrino energy.) 

\ The recent Kamiokandc upper limit of 4 x 1CP 14 muons cm" 2 s _1 sr _1 for zenith angles larger 
than 150° [17] is slightly less stringent than the 1MB limit we consider, which corresponds 
to 3.7 x 1CT 14 muons cm" 2 s" 1 sr" 1 . 




(4.1) 
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For a simplified model of the Earth with uniform density p© = 5.5 gem -3 and 
radius i?© = 6.37 x 10 8 cm, the effective aperture is (neglecting the depth of the 
underground detector relative to -R©), 

ni(E v .) = J dQexp[- 2^(^)1 cos 0(-cos0), (4.2) 

where the integral extends over the geometrical aperture of the detector, d is the 
zenith angle and ki{E Vi ) is the neutrino absorption coefficient in the Earth, given 
by 

h{E Vi ) = £^o Vi K{E Vi ), (4.3) 
m N 

with mN the nucleon mass and a Vi ^{E Vi ) the total neutrino-nucleon cross section. 
For the neutrino energies under consideration, <j v j$(E Vi ) includes only the charged 
current cross section <j^^(E Vi ), because the energy and momentum fractions trans- 
ferred to the nucleon in a neutral current process are negligible at these energies. 

The charged current cross section a^(E Ui ) is well-known for E Vi ^ lOTeV: 

^ = 0.67xl0- 35 cm 2 (^) , (4.4) 



for a neutrino and 



^ = 0.34 x 10~ 3 W , (4.5) 



for an antineutrino [5]. At higher energies the charged current cross section becomes 
more and more uncertain — by as much as a factor of 10 at E Vi ~ 10 9 TeV - 
because of the poor knowledge of nucleon structure functions at small arguments 
[11]. For this reason, we have not attempted a precise calculation of a^^(E Ui ) 
from a set of theoretical structure functions. We have used the differential charged 
current cross section up to E Vi = 10 7 TeV given in ref. [8]. At still higher energies, 
we have matched the asymptotic form of the cross section in ref. [7] to the results 
of ref. [8]. 
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In figure 3 we show the effect of absorption in the Earth by plotting the effective 
aperture (4.2) integrated below the horizon: 



nr™(E„j 2:7,7 



1 — exp 



(4.6) 



with a© = mN/2i? e p e = 2.4 x 1CT 34 cm 2 . This effective aperture differs very little 
from the one obtained in ref. [19] using a more elaborate model of the Earth. As 
we see from the figure, the Earth severely attenuates the flux of neutrinos of energy 
exceeding ~ 10 5 TeV, becoming nearly opaque at ~ 10 10 TeV. 

Note that for 10~ 33 cm 2 , i.e. at E Vi ^ 10 3 TeV, the effective aperture from 
below the horizon is inversely proportional to the neutrino-nucleon scattering cross 
section, 

^ el ° w (^) ^ 2tt J . (4.7) 

The resonant reaction u e e~ — > W~ — > "anything" severely depletes the u e 
flux from below the horizon at energies around E„ c ~ 7 x 10 3 TeV [20]. However 
since we do not assume any predominant neutrino type in the decay neutrino flux, 
the u e flux from below accounts for only one eighth of the total rate of contained 
events. It is therefore a reasonable approximation, for our purposes, to neglect this 
resonance. 

We present now the transfer functions Pi(E Vi ) for the experimental data sets 
under consideration. For contained events we have 

mm(E cut ,E Vi ) ^ 

PT A {E Vi ) = 0{E Vi -E^)^ J dE h ^, (4.8) 

E th 

where E t ^ is the experimental energy threshold for the lepton energy E[., E cu t 
is an experimental cutoff (2.5 GeV for 1MB and infinite for Frejus), M is the 
detector mass and A^ nuc i = 6.02 x 10 32 (M/kton) is the number of nucleons in 
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the detector. If the charged current cross section is written in units of 10 38 cm 2 , 
a^(E Vi ) = (7^38 x lCT 38 cm 2 , then the transfer function for contained events is 

nan(E ca t,E„ i ) 

Pr\E Vi ) ~ 6.0 x 10- 6 cm 2 kton- 1 9{E Vi - E th ) f dE k ^^. (4.9) 

E th 

Both Frejus and 1MB data sets include neutrinos coming from all solid angles and 
their effective aperture computed from eq. (4.2) varies from An to 2n as the energy 
is increased; the neutrino flux is reduced at most by a factor of 2 at the highest 
energies. 

In 1MB, and other water-Cerenkov detectors, there is also the possibility that 
the hadronic fragments produced in the neutrino-nucleus collision give a detectable 
amount of Cerenkov light. Their contribution ^/hadr(^) snom d then be added 
to eq. (4.9). In the appendix, we present an estimate of the contribution from 
such 'hadronic blasts' and show that this is important only for very energetic neu- 
trinos, -E^^10 7 TeV, where, however, the signal from EAS's gives more stringent 
constraints. 

The product Pf ont (E Vi ) fi? ont (£' I/i ) thus obtained for the Frejus and 1MB con- 
tained events is shown in figure 4(a) for neutrinos (solid lines) and antineutri- 
nos (dotted lines) as function of the neutrino (or antineutrino) energy E Vi . The 
units are chosen such that the vertical axis directly gives the number of events 
per kiloton-year corresponding to a unit neutrino flux of 1 cm _2 s _1 sr _1 . For 
E v . k, 10 6 TeV we calculate Pf ont (E Vi ) Vq orA (E Vi ) ~ 3.8 x 10~ 5 cm 2 sr kton" 1 a i)38 
for the Frejus detector. The 1MB curve (curve 1) above 10 7 TeV is due to 'hadronic 
blasts' as discussed in the appendix. 

The transfer function for the flux of up-going muons is (see ref. [21]) 

oo E Vi oo 

P?{E Vi ) = Q{E V% - E th ) J dE^ J dE, J dXg(X, E'^ E,) (4.10) 

E th E'„ 

for V{ = v^v^, and Pf{E Vi ) = for the other neutrino types. Here E^ and 
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E'p are the muon energies at production and at the detector respectively, and 
X = /p ro ck/ m N is the column density of rock, i.e. the number of nucleons per 
unit area encountered by a muon travelling a length / in rock. Notice that Pf is 
adimensional. The probability that a muon with initial energy has an energy 
between E'^ and E'^ + dE'^ after traversing an amount x of rock is denoted by 
g (X, E'^, E^) dE'^. We assume a uniform rock density of p roc k = 2.6 gem -3 in the 
region surrounding the detector. Following ref. [6], we make the approximation 
that the final muon energy E'^ coincides with its mean value (with no dispersion): 

El = (Bp + e)e~^ x - e, (4.11) 

with e ~ 0.51 TeV and 7 _1 = 1.54 x 10 29 cm~ 2 . The integral over E'^ in eq. (4.10) 
can then be performed, and we obtain 

E v . 

P?{E n ) = 6{E Vi - E th ) J dE,X th (E,)^, (4.12) 
for v i = v il ,v il . Here 

w^'Kot) (413) 

is the column density traversed by muons produced with energy E^ which reach 
the detector with threshold energy E t ^. A plot of Pjf{E Vi ) times Q?(E Vi ), obtained 
by integrating eq. (4.2) over zenith angles larger than 98°, is shown in figure 4(b). 
The decrease of P^Vt^ for E Vi H 10 7 TeV is due to absorption by the Earth. 

The final signal we consider is the rate of EAS's per unit solid angle. Its 
transfer function is 

P f KS ( E n) = ^J {Ev f uMEuME Vi - E th ), (4.14) 

where the index i stands for u e and u T , which can generate electromagnetic or 
hadronic showers in the atmosphere. Muons from charged current z/ M -nucleon in- 
teractions do not trigger electromagnetic showers, since their radiation length for 
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bremstrahlung in air (10 5 g/cm 2 ) is much larger than the atmosphere thickness 
(1030 g/cm 2 ). The product Pf ,AB {E Vi ) Qf AS (E Vi ) is shown in figure 4(c) for the 
four experimental energy thresholds of the Fly's Eye detector [18]. 

We are now in a position to compare the expected signals to the experimental 
limits. 

5. Present constraints 

We assume that the same numbers of neutrinos and antineutrinos of each type, 
B Vc = B Uc = B v = B Vti = ... = B u , are produced in x decays and that their 
production energy is always E e = ^m x . Inserting the decay-generated neutrino 
flux, eq. (3.5), into eq. (4.1) and using the appropriate transfer functions and 
effective apertures described in section 3, we obtain the expected signals in terms 
of the decay lifetime r x and the quantity B ll m x Y Xp = B l/ m x Y Xo e t °/ Tx which is 
proportional to the primordial energy density of the decaying particles. Notice 
that for r x ^> to, this quantity is ~ 25.5 eV [B v f2 XQ h 2 ), where Q Xo is the present x 
mass density in units of the critical cosmological density. 

We present the results in figure 5. The shaded regions are excluded by the 
present experimental data. The solid lines refer to the limit on upward-going muons 
from 1MB, the dotted and short-dashed lines to the Frejus and 1MB contained 
events, respectively, and the long-dashed lines to the Fly's Eye EAS's. These 
are essentially bounds on the relic energy density of the decaying particle, taking 
B v — 1; for B v < 1, these lines are to be proportionally shifted upward. As noted 
earlier, the experimental limits on contained events can, in principle, be improved 
by a factor of ~ 10 if the signal due to atmospheric neutrinos is accounted for; the 
corresponding bounds should then be scaled downwards by the same factor. The 
short-dashed-dotted line corresponds to a present mass density QqIi 2 = 1 either in 
x particles (for r x ^to) or m its decay products (for r x &to); the region Q^h? > 1 
is excluded by the observational lower limits to the age and present expansion rate 
of the universe. For comparison we also show as a long-dashed-dotted line (in the 
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lower left quadrant) the upper bound on m x Y Xp for unstable particles decaying 
into electromagnetically interacting particles. This is obtained by requiring that 
the abundance of the primordially synthesised light elements D, 3 ife, ^He and 
7 Li not be excessively altered from their observationally inferred values by the 
electromagnetic cascades initiated by the decay products [22]. In fact this bound 
also applies to unstable particles decaying into neutrinos, since the decay neutrinos 
can initiate similar electromagnetic cascades through the process vv — > e + e~, 
where the target (anti)neutrinos belong to the thermal background. (This has 
also been considered in ref. [23]; however these authors do not calculate cascade 
generation correctly and obtain an overly restrictive bound.) 

Figures 5(a-d) correspond to m x = 1, 10 5 , 10 6 and 10 10 TeV respectively. 
As the x mass increases, the bounds at r x & to first shift to the left as the decay 
neutrinos become more energetic and the signals go further above the experimental 
thresholds. Then they proceed to move to the right because the neutrino absorption 
redshift decreases with increasing neutrino energy, hence the neutrinos produced 
by relatively short-lived particles do not survive until the present. The mass at the 
turning point is given approximately by solving m x ~ 3-Eth[l + z a (^m x )\ with the 
help of eqs. (2.21); its value is 1 x 10 3 TeV and 1 x 10 2 TeV for contained events in 
Frejus and 1MB respectively, 8 x 10 2 TeV for the 1MB upward-going muons, and 
1 x 10 8 TeV, 1 x 10 9 TeV, 5 x 10 9 TeV and 3 x 10 10 TeV for the four Fly's Eye 
thresholds. 

When r x ^ to, the best bounds on the relic abundance of the decaying particles 
come from the 1MB limit on upward-going muons at m x ^5 x 10 5 TeV and the 
Fly's Eye limits on EASs at m x ^5 x 10 5 TeV. We can invert the argument and 
consider the interesting case f2 Xo h 2 ~ 1, i.e. when the x particles are assumed to 
constitute the dark matter today* A corresponding lower bound on its lifetime 

* In this case the actual spatial distribution of the relic particles should be taken into account, 
e.g. their likely concentration in the halo of our Galaxy. This would yield even stricter 
constraints. Preliminary work has been reported in ref. [24] and a more detailed study is in 
progress. 
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versus its mass can then be inferred and is plotted in figure 6. For m x ^ 30 TeV, 
this bound gets stronger with increasing m x as the mean energy of the decay 
generated neutrinos rises over the 1MB detection threshold. For m x ^ 30 TeV, the 
lower bound on t x is inversely proportional to m x (cf. eq. (3.9)), apart from the 
jumps at ~ 10 5 TeV, ~ 10 6 TeV, ~ 10 7 TeV and ~ 10 8 TeV corresponding to the 
Fly's Eye energy thresholds. No bound exists for m x ^5 x 10 14 TeV, since the 
universe is opaque to such high energy neutrinos at the present epoch. 

6. Conclusions 

We have considered constraints on the lifetime, the abundance and the mass of 
unstable relic particles decaying into neutrinos at cosmological epochs, taking into 
account that both the early universe and the Earth are opaque to very high energy 
neutrinos. We have evaluated the signals expected from a diffuse background 
of decay-generated neutrinos in underground nucleon decay experiments and at 
cosmic ray observatories. Comparing these to the present limits on the flux of 
non-atmospheric neutrinos, we find severe bounds on the relic abundance of such 
heavy particles; in particular, such particles must be very long-lived indeed in order 
to constitute the dark matter today. These bounds are of relevance to massive 
metastable particles such as technicolour baryons and 'cryptons' (bound states in 
the hidden sector of superstring-inspired models) as discussed elsewhere [22] . 
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Recently, we became aware of ref. [25], where the detection of neutrinos from 
cosmic relic particles is also studied, in particular the effects due to absorption in 
the early Universe. However this work assumes that the decaying particles were 
thermally produced (with a calculable abundance) in the early Universe, whereas 
we have presented results in a general form, applicable to any relic particle. This, 
in fact, is essential in order to consider particles with masses over a few hundred 
TeV. We also believe that we have addressed experimental issues in more detail. 



APPENDIX 

Here we present an estimate of the contribution of 'hadronic blasts' to the 1MB 
transfer function for contained events. We calculate the visible energy equivalent 
to the Cerenkov light output from such a 'blast' and then compare it to the visible 
energy observable in the 1MB detector. 

Let W be the energy transferred to the nucleus in the neutrino interaction. This 
is also presumably the energy available to the hadronic shower generated by the 
nuclear fragments. The visible energy E v i s is defined to be the energy of a fictitious 
initial electron generating an electromagnetic shower with the same Cerenkov light 
output [13]. The physics of electromagnetic showers [26] then relates E v i s to the 
'detectable' track length X^, 

where E c is the critical energy separating the domains where ionization and radi- 
ation energy losses dominate and is approximately given by 

SOOMeV 

Z+1.2 K J 

Here Xq(— 36.08 g/cm 2 [27]) is the electron radiation length in water and F(z) is 
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approximately 

F(z)^e*(l + zln_L_), (A.3) 

with 

, r- Q Z Ed Ed 

z ~ 4.58 ~ A. 4 

A E c 28MeV V ; 

(the numerical values refer to water). In the 1MB detector, Ed = 1.52 m e + 
30MeV + 140 MeV = 170 MeV for an electron [14], hence eq. (A.l) reads 

E vis = 4.8 x 10^ 4 MeV X d (W) (A.5) 

with Xd(W) in g/cm 2 . 

It remains now to estimate Xd(W), the path length of charged particles in the 
hadronic shower with energy above the detection threshold. The mean number of 
charged particles n^W) as a function of the available energy W has been studied 
in deep inelastic scattering and incorporated into the Lund Monte Carlo program 
[28]. By fitting fig. 8 of this reference we obtain: 



n<&(W) = 1.67 + 0.211 exp 



(A.6) 



We assume now that all charged particles in the first generation of the shower are 
above the Cerenkov threshold. This is a good approximation at the high energies 
that will turn out to be relevant for hadronic blasts. With this assumption, Xd(W) 
is the product of n^iW) and the mean path length of a charged particle. A lower 
bound to the latter is one nuclear interaction length, X nuc i = 84.9 g/cm 2 in water. 
This leads to a lower bound for Xd(W): 

X d (W) > X min (W) ~ n ch (W)X nucl . (A.7) 

An upper bound is obtained by multiplying X m [ n (W) by the typical length of a 
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hadronic shower in units of interaction lengths 

W 



X d (W) < X max (W) ~ n ch (W)X nucl 



5.45 + 0.89 In 



GcV 



(A.i 



A lower and an upper bound to the visible energy E Y - m (W) can then be obtained 
from eq. (A. 5). 

This range then has to be compared with the 1MB energy threshold for 7r ± 
detection, E n = 1.52 + 30MeV + 140 MeV = 382 MeV, and with the highest 
energy analyzed, E Y - m = 2500 MeV [13]. The result of such a comparison is that for 
W^IO 2 TeV there is probably not enough Cerenkov light for the 'blast' to be seen, 
while for W^10 6 TeV the detector is probably overloaded (more than 900 PMTs 
fired [13]) or otherwise not sufficiently efficient to detect the blast. However for 
10 2 TeV^H/^10 6 TeV, such blasts, if they do occur, should already be present in 
the sample of ref. [13], probably as multiple-ring events. 

However, the neutrino energy E v required to have W^IO 2 TeV is £^10 7 TeV. 
This comes from the kinematic relation W 2 ~ 2m p E l/ (l — x)y, with m p the nucleon 
mass and x and y the usual deep inelastic scattering variables, together with the 
consideration that at high energies the W or Z propagator restricts the important 
values of x and y to x ~ and y — 1. So 'hadronic blasts' in Cerenkov detectors 
turn out to be important only for very high energy neutrinos, where (at least for the 
purposes of this paper) there already are much better bounds on cosmic neutrino 
fluxes from EAS arrays. 

For the sake of completeness, we show in fig. 4a (curve 1 at high energies) the 
contribution of hadronic blasts in 1MB, obtained by numerical integration of: 



2m p E v WL* 
Annuel f , ^.9 f i ti^2 ' 



^(^) = ^ / dQ 2 J dW 



w 2 . 

mm 



(2m p E v y y 



da CC da NC- 

7nr + 7nr ■ A - 9 

axay axay 



Here Q 2 ~ 2m p E u xy, W m [ n = 10 2 TeV, W max = 10 6 TeV and the contributions from 
charged and neutral currents have been summed. The EHLQ structure functions 
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have been used together with a McKay-Ralston asymptotic form at low x (as in 
ref. [10]). 
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FIGURE CAPTIONS 



Fig. 1. The absorption redshift z K (line 1) for cosmic neutrinos as a function of the 
neutrino energy at emission E e taking QqIi 2 = 1. The other lines indicate: 
(2) the boundary between the regions where absorption due to annihilation 
and scattering dominate; (3) the present epoch; (4) the Z boson pole; (5) the 
epoch of matter-radiation equality; (6) the epoch of light neutrino decoupling. 

Fig. 2. The present energy spectrum of decay generated neutrinos for r x = 10~ 5 to 
and 3E e = 10 5 TeV (line 1), 10 7 TeV (line 2) and 10 9 TeV (line 3). The full 
lines show the effects of cosmological neutrino absorption. 

Fig. 3. The effective detector aperture, integrated below the horizon, for neutrinos 
(solid line) and antineutrinos (dotted line), demonstrating the opacity of the 
Earth at high energies. 

Fig. 4. The product of the transfer function Pi(E Ui ) and of the effective aperture 
Qi(E Ui ) for the three experimental data sets we consider: (a) contained events 
in 1MB (curve 1) and Frejus (curve 2) (b) 1MB upward-going muons and (c) 
Fly's Eye EAS's (four thresholds). The dotted lines corresponds to antineu- 
trinos. 

Fig. 5. The bound on the primordial energy density of the decaying particle mul- 
tiplied by the branching ratio into neutrinos, B iy m x Y Xp , as a function of its 
lifetime t x , for various choices of its mass m x . The shaded regions are ex- 
cluded by the present experimental data. The various lines refer to: 1MB 
upward-going muons (solid lines), Frejus and 1MB contained events (dotted 
and short-dashed lines respectively), Fly's Eye EAS's (long-dashed lines). 
Also indicated are the upper bound on the total energy density QqIi 2 = 1 
(short-dashed-dotted line) and the upper bound inferred from considerations 
of primordial light element abundances (long-dashed-dotted line) . 

Fig. 6. The lower bound on the x particle lifetime versus its mass for a present 
density (2 Xo h 2 = 1 in the relic particles, assuming unit branching ratio into 
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